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New optical arrangements with two single-mode input fibers and a fiber-optic coupler are devised 
to measure the instantaneous velocity difference and local velocity. The fibers and the coupler are 
polarization-preserving to guarantee a high signal-to-noise ratio. When the two input fibers are 
used to collect the scattered light with the same momentum transfer vector but from two spatially 
separated regions in a flow, the obtained signals interfere when combined via the fiber-optic coupler. 
The resultant light received by a photomultiplier tube contains a cross-beat frequency proportional to 
the velocity difference between the two measuring points. If the two input fibers are used to collect 
the scattered light from a common scattering region but with two different momentum transfer 
vectors, the resultant light then contains a self-beat frequency proportional to the local velocity at 
the measuring point. The experiment shows that both the cross-beat and self-beat signals are large 
and the standard laser Doppler signal processor can be used to measure the velocity difference and 
local velocity in real time. The new technique will have various applications in the general area of 
fluid dynamics. 

OCIS codes: 120.7250, 060.2420, 060.1810, 030.7060. 



I. INTRODUCTION 



Measurements of the velocity difference or the relative 
velocity, 5v(£) = v(x + £) — v(x), between two spatial 
points separated by a distance £ have important appli- 
cations in fluid dynamics. For example, in the study of 
turbulent flows one is interested in the scaling behavior 
of 8v(£) over varying distance £, when £ is in the in- 
ertial range, in which the kinetic energy cascades at a 
constant rate without dissipation. [fj],p| If the separation 
I is smaller than the Kolmogorov dissipation length, jl],^| 
the measured 5v(£) becomes proportional to the velocity 
gradient dv/dr ~ 5v(£)/£ (assuming £ is known), a useful 
quantity which is needed to determine the energy dissi- 
pation and flow vorticity. In many experimental studies 
of fluid turbulence, one measures the local velocity as 
a function of time and then uses Taylor's frozen turbu- 
lence assumption to convert the measured temporal vari- 
ations into the spatial fluctuations of the velocity field. 
H The frozen turbulence assumption is valid only when 
the mean velocity becomes much larger than the velocity 
fluctuations. For isotropic turbulent flows with a small 
mean velocity, direct measurement of Sv(£) is needed. 



Over the past several years, the present authors and 
their collaborators have exploited the technique of ho- 
modyne photon correlation spectroscopy (HPCS) to mea- 
sure 5v(£). |||,5j With the HPCS scheme, small particles 
seeded in a flowing fluid are used to scatter the inci- 
dent laser light. The scattered light intensity I(t), which 
fluctuates because of the motion of the seed particles, 
contains Doppler beat frequencies of all particle pairs in 
the scattering volume. For each particle pair separated 
by a distance £ (along the beam propagation direction), 
their beat frequency is Au>2 = q • Sv(£), where q is the 
momentum transfer vector. The magnitude of q is given 
by q = (47m/A) sin(0/2), where 9 is the scattering angle, 
n is the refractive index of the fluid, and A is the wave- 
length of the incident light. Experimentally, the Doppler 
beat frequency Au>2 is measured by the intensity auto- 
correlation function, pj 



9(r) 



(I(t + T)I(t)) 
(I(t)) 2 



1 + 6G(r), 



(1) 



where 6 (< 1) is an instrumental constant and henceforth 
we set 6=1. The angle brackets represent a time average 
over t. 
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It has been shown that G(t) in Eq. ([!]) has the form 



G(r) = / dr h(r) dSv P(Sv, r) cos{qSvr), (2) 

JO J -oo 

where Sv is the component of 8\ in the direction of 
q, P(Sv,r) is the probability density function (PDF) of 
6v(r), and h(r)dr is the number fraction of particle pairs 
with separation r in the scattering volume. Equation 
(||) states that the light scattered by each pair of par- 
ticles contributes a phase factor cos(qr5v) (because of 
the Doppler beat) to the correlation function G(r), and 
G(t) is an incoherent sum of these ensemble averaged 
phase factors over all the particle pairs in the scattering 
volume. In many previous experiments, [fz| — |l0|| the length 
L of the scattering volume viewed by a photodetector was 
controlled by the width S of a slit in the collecting optics. 

While it is indeed a powerful tool for the study of tur- 
bulent flows, the HPCS technique has two limitations 
in its collecting optics and signal processing. First, a 
weighted average over r is required for G(r) because the 
photodetector receives light from particle pairs having a 
range of separations (0 < r < L). As a result, the mea- 
sured G(r) contains information about <5v(£) over various 
length scales up to L. With the single slit arrangement, 
the range of L which can be varied in the experiment is 
limited. The lower cut-off for L is controlled by the laser 
beam radius a. The upper cut-off for L is determined by 
the coherence distance (or coherence area) at the detect- 
ing surface of the photo-detector, over which the scat- 
tered electric fields are strongly correlated in space. Q 
When the slit width S becomes too large, the photode- 
tector sees many temporally fluctuating speckles (or co- 
herence areas) , and consequently fluctuations in the scat- 
tered intensity I(t) will be averaged out over the range 
of g-values spanned by the detecting area. 

Recently, we made a new optical arrangement for 
HPCS, with which the weighted average over r in Eq. 
(^|) is no longer needed and the upper limit for L can be 
extended to the coherence length of the laser. In the ex- 
periment, S two single mode, polarization-maintaining 
(PM) fibers are used to collect light with the same po- 
larization and momentum transfer vector q but from two 
spatially separated regions in a flow. These regions are 
illuminated by a single coherent laser beam, so that the 
collected signals interfere when combined using a fiber- 
optic coupler, before being directed to a photodetector. 
With this arrangement, the measured G(r) becomes pro- 
portional to the Fourier cosine transform of the PDF 
P(Sv,r). 

The second limitation of HPCS is related to signal pro- 
cessing. The correlation method is very effective in pick- 
ing up small fluctuating signals, but the resulting corre- 
lation function G(r) is a time-averaged quantity. There- 
fore, the correlation method is not applicable to unstable 
flows. Furthermore, information about the odd moments 
of P(Sv, r) is lost, because the measured G(r) is a Fourier 



cosine transform of P(6v,r). 

In this paper, we present a further improvement for 
HPCS, which is free of the two limitations discussed 
above. By combining the new fiber-optic method with 
the laser Doppler velocimetry (LDV) electronics, we are 
able to measure the instantaneous velocity difference 
5v(£, t) and local velocity v(x, t) at a high sampling rate. 
With this technique, the statistics of 5v(£, t) and v(x, t) 
are obtained directly from the time series data. The 
new method of measuring v(x, t) offers several advan- 
tages over the standard LDV. The remainder of the paper 
is organized as follows. In Section 2 we describe the ex- 
perimental methods and setup. Experimental results are 
presented and analyzed in Section 3. Finally, the work is 
summarized in Section 4. 



II. EXPERIMENTAL METHODS 

A. Measurement of the velocity difference 

Figure 1 shows the optical arrangement and the flow 
cells used in the experiment. A similar setup has been 
described elsewhere, j^] and here we mention only some 
key points. 




FIG. 1. (a) Scattering geometry for the velocity difference 
measurement, ki, incident wave vector; k s , scattered wave 
vector; and q = k s — ki. (b) Experimental setup for the ve- 
locity difference measurement in rigid body rotation, (c) Flow 
cell and optical arrangement for a jet flow. 

As shown in Fig. 1(b), an incident beam from a 
Nd:YV04 laser with a power range of 0.5-2W and a wave- 
length of A=532 nm is directed to a flow cell by a lens. 
With the aid of a large beam splitting cube, two single- 
mode, polarization-maintaining (PM) fibers collect the 
scattered light from two different spots along the laser 
beam with a separation L. The two PM fibers are con- 
nected to a fiber-optic coupler (purchased from OZ Op- 
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tics fLl| ), which combines the light from the two input 
fibers and split the resultant light evenly into two output 
fibers. A graded index lens is placed at each end of the 
fiber to collimate the light entering (or exiting from) the 
fiber core. Each input fiber is mounted on a micrometer- 
controlled translation stage and the distance L can be 
adjusted in a range of 0-25 mm in steps of 0.01 mm. The 
output fibers of the optical coupler are connected to two 
photomultiplier tubes (PMT1 and PMT2). PMT1 is op- 
erated in the digital mode and its output signal is fed to a 
digital correlator (ALV-5000). PMT2 is operated in the 
analogue mode and its output signal is fed to a LDV sig- 
nal processor (TSI IFA655). An oscilloscope connected 
to PMT2 directly views the analogue signals. A low- noise 
preamplifier (Stanford Research SR560) further amplifies 
the analogue output of PMT2 before it goes to the LDV 
signal processor. 

As shown in Fig. 1(a), the electric fields detected by 
each input fiber sum in the coupler and consequently in- 
terfere. In the experiment, we obtain the beat frequency, 
Alu 2 — q v! q v 2 = q-Sv(L), in two different ways. One 
way is to measure the intensity auto-correlation function 
g(r) in Eq. (|l]). With the ALV correlator, it takes ~ 1 
minute to collect the data with an adequate signal-to- 
noise ratio. The other way is to use the LDV signal pro- 
cessor to measure the instantaneous beat frequency Aw2, 
giving velocity differences in real time. The LDV signal 
processor is essentially a very fast correlator and thus re- 
quires the beat signal to be large enough so that no signal 
averaging is needed. In the experiment to be discussed 
below, we use both methods to analyze the beat signals 
and compare the results. 

It has been shown that the correlation function g(r) 
has the form: Q 

■' M = 1+ (#W G ' (r) + (TTW ft(T) 

= 1 + 6,G,(t)+6«G c (t), (3) 

where I\ and I2 are the light intensities from the two 
input fibers. When 1\ = I2, one finds b s = b c = 0.5. 
If one of the input fibers is blocked (i.e., I2 — 0), we 
have g{r) = 1 + G s (r), where G s (t) is the self-beat cor- 
relation function for a single fiber. When the separation 
L between the two input fibers is much larger than the 
spot size viewed by each fiber, the cross-beat correlation 
function G c (r) takes the form 

/+00 
dSv P{6v)cos(q5v(L)T). (4) 
-00 

Two flow cells are used in the experiment. The first 
one is a cylindrical cuvette having an inner diameter of 
2.45 cm and a height of 5 cm. The cuvette is top mounted 
on a geared motor, which produces smooth rotation with 
an angular velocity u) = 2.5 rad/s. The cell is filled with 
1,2-propylene glycol, whose viscosity is 40 times larger 
than that of water. The whole cell is immersed in a large 



square index-matching vat, which is also filled with 1,2- 
propylene glycol. The flow field inside the cell is a simple 
rigid body rotation. With the scattering geometry shown 
in Fig. 1(b), the beat frequency is given by Aw2 = k s u>L 
with k s = 27m/ A. The sample cell is seeded with a small 
amount of polystyrene latex spheres. For the correla- 
tion measurements, we use small seed particles of 1.0 [xm 
in diameter. By using the small seed particles, one can 
have more particles in the scattering volume even at low 
seeding densities. This will reduce the amplitude of the 
number fluctuations caused by a change in the number of 
particles in each scattering volume. The particle number 
fluctuations can produce incoherent amplitude fluctua- 
tions to the scattered light and thus introduce an extra 
(additive) decay to g(r). [|l6) Large particles 4.75 fim in 
diameter are used to produce higher scattering intensity 
for instantaneous Doppler burst detection. Because the 
densities of the latex particles and the fluid are closely 
matched, the particles follow the local flow well and they 
do not settle much. 

The second flow cell shown in Fig. 1(c) is used to gen- 
erate a jet flow in a 9cm x 9cm square vat filled with 
quiescent water. The circular nozzle has an outlet 2 mm 
in diameter and the tube diameter before the contraction 
is 5.5 mm. The nozzle is immersed in the square vat and 
a small pump is used to generate a jet flow at a constant 
flow rate 0.39 cm 3 /s. The mean velocity at the nozzle 
exit is 12.4 cm/s and the corresponding Reynolds num- 
ber is Re=248. This value of Re is approximately three 
times larger than the turbulent transition Reynolds num- 
ber, Re c ~ 80, for a round jet. jl2| Because of the large 
area contraction (7.6:1), the bulk part of the velocity pro- 
file at the nozzle exit is flat. This uniform velocity pro- 
file vanishes quickly, and a Gaussian-like velocity profile 
is developed in the downstream region, 2-20 diameters 
away from the nozzle exit. |D| As shown in Fig. 1(c), 
the direction of the momentum transfer vector q (and 
hence the measured velocity difference) is parallel to the 
jet flow direction, but the separation L is at an angle of 
45° to that direction. 



B. Measurement of the local velocity 

In the measurement of the velocity difference, we use 
two input fibers to collect the scattered light with the 
same q (i.e. at the same scattering angle) but from two 
spatially separated regions in the flow. We now show 
that with a different optical arrangement, the fiber-optic 
method can also be used to measure the instantaneous 
local velocity v(x, t). Instead of collecting light from two 
spatially separated regions with the same q, we use the 
two input fibers to collect light from a common scatter- 
ing volume in the flow but with two different momentum 
transfer vectors qi and q2 (i.e. at two different scatter- 
ing angles). Figure 2(a) shows the schematic diagram 
of the scattering geometry. The collected signals at the 
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two scattering angles are combined by a fiber-optic cou- 
pler, and the resultant light is modulated at the Doppler 
beat frequency: Q Ami = qi • v q 2 • v = Aq ■ v(x), 
where Aq = qi — q 2 . The magnitude of Aq is given 
by Aq = (Ann/X) sin(a/2), with a being the acceptance 
angle between the two fibers. 

The principle of using the scattered light at two dif- 
ferent scattering angles to measure the local velocity has 
been demonstrated many years ago. 14 What is new 
here is the use of the fiber-optic coupler for optical mix- 
ing. The fiber-optic technique simplifies the standard 
LDV optics considerably. As shown in Fig. 2(b), in the 
standard LDV arrangement the two incident laser beams 
form interference fringes at the focal point. When a seed 
particle traverses the focal region, the light scattered by 
the particle is modulated by the interference fringes with 
a frequency, jl5) Aui = qi • v q 2 ■ v = Aq • v(x). 
The magnitude of Aq has the same expression Aq = 
(47rn/A) sin(a/2) as shown in the above, but a now be- 
comes the angle between the two incident laser beams. 
The main difference between the standard LDV and the 
new fiber-optic method is that the former employs two 
incident laser beams and a receiving fiber [Fig. 2(b)], 
while the latter uses only one incident laser beam and 
two optical fibers to measure each velocity component 
[Fig. 2(a)]. Consequently, the beat frequency Alui in 
Fig. 2(a) is independent of the direction (ki) of the inci- 
dent laser beam, whereas in Fig. 2(b) it is independent 
of the direction (k s ) of the receiving fiber. 
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FIG. 2. (a) One-beam scattering geometry for the local ve- 
locity measurement, ki, incident wave vector; (k s )i and (k s )2, 
two scattered wave vectors; qi = (k s )i — ki; q2 = (k s )2 — ki; 
Aq = q2 — qi. (b) Two-beam scattering geometry for the 
local velocity measurement. (ki)i and (ki)2, two incident 
wave vectors; k s , scattered wave vector; qi = k s — (ki)i; 
q2 = k s — (ki)2; Aq = q2 — qi. (c) Schematic diagram of a 
one-beam probe for the local velocity measurement. S, mea- 
suring point; P, lens; M, frequency modulator; C, fiber-optic 
coupler. 

With the one-beam scheme, one can design various op- 
tical probes for the local velocity measurement. Figure 



2(c) shows an example, which would replace the commer- 
cial LDV probe by reversing the roles of the transmitter 
and receiver. The two input fibers aim at the same mea- 
suring point S through a lens P, which also collects the 
back-scattered light from S. The frequency modulator M 
shifts the frequency of the light collected by the input 
fiber 1, before it is combined via the fiber-optic coupler 
C with the light collected by the input fiber 2. The re- 
sultant light from an output fiber of the coupler contains 
the beat frequency Awi and is fed to a photodetector. 
Because the measured Au>\ is always a positive number, 
one cannot tell the sign of the local velocity when zero 
velocity corresponds to a zero beat frequency. The fre- 
quency shift by the modulator M causes the interference 
fringes to move in one direction (normal to the fringes) 
and thus introduces an extra shift frequency to the mea- 
sured beat frequency. This allows us to measure very 
small velocities and to determine the sign of the mea- 
sured local velocity relative to the direction of the fringe 
motion (which is known). Jl5[ The other output fiber of 
the coupler can be used as an alignment fiber, when it is 
connected to a small He-Ne laser. With the reversed He- 
Ne light coming out of the input fibers, one can directly 
observe the scattering volume viewed by each input fiber 
and align the fibers in such a way that only the scattered 
light from the same measuring point S is collected. 

The one-beam probe has several advantages over the 
usual two-beam probes. To measure two orthogonal ve- 
locity components in the plane perpendicular to the in- 
cident laser beam, one only needs to add an extra pair of 
optical fibers and a coupler and arrange them in the plane 
perpendicular to that shown in Fig. 2(c) (i.e., rotate the 
two-fiber plane shown in Fig. 2(c) by 90°). With this 
arrangement, the four input fibers collect the scattered 
light from the same scattering volume but in two orthog- 
onal scattering planes. Because only one laser beam is 
needed for optical mixing, a small single- frequency diode 
laser, rather than a large argon ion laser, is sufficient for 
the coherent light source. In addition, the one-beam ar- 
rangement does not need any optics for color and beam 
separations and thus can reduce the manufacturing cost 
considerably. With the one-beam scheme, one can make 
small invasive or non-invasive probes consisting of only 
three thin fibers. One can also make a self-sustained 
probe containing all necessary optical fibers, couplers, 
photodetectors, and a signle-frequency diode laser. 



III. RESULTS AND DISCUSSION 

A. Velocity difference measurements 

We first discuss the measurements of the velocity dif- 
ference in rigid body rotation. Figure 3(a) shows a sam- 
ple oscilloscope trace of the analogue output from PMT2 
when the separation L=1.0 mm. The signal is amplified 
1250 times and band-pass-filtered with a frequency range 
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of 1-10 kHz. This oscilloscope trace strongly resembles 
the burst signals in the standard LDV. The only differ- 
ence is that the signal shown in Fig. 3(a) results from the 
beating of two moving particles separated by a distance 
L. Figure 3(a) thus demonstrates that the beat signal 
between the two moving particles is large enough that a 
standard LDV signal processor can be used to measure 
the instantaneous velocity difference 5v(£, t) in real time. 

Figure 3(b) shows the measured beat frequency Aw2 
as a function of separation L. The circles are obtained 
from the oscilloscope trace and the triangles are obtained 
from the intensity correlation function g(r). The two 
measurements agree well with each other. The solid 
line is the linear fit Aw2 = 41.78L (10 3 rad/s), which 
is in good agreement with the theoretical calculation 
Aw 2 = k s ujL = 42.28L (10 3 rad/s). This result also 
agrees with the previous measurements by Du et al. || 
Because Aw2 increases with L, one needs to increase the 
laser intensity at large values of L in order to resolve A0J2 ■ 
The average photon count rate should be at least twice 
the measured beat frequency. 
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FIG. 3. (a) Oscilloscope trace of a typical beat burst be- 
tween two moving particles separated by a distance L=1.0 
mm. The signal is obtained in rigid body rotation, (b) Mea- 
sured beat frequency Au>2 as a function of separation L. The 
circles are obtained from the oscilloscope trace and the trian- 
gles are obtained from the intensity correlation function g(r). 
The solid line shows a linear fit to the data points. 

We now discuss the time series measurements of 
Sv(L, t) in a jet flow using the LDV signal processor. The 
jet flow has significant velocity fluctuations as compared 
with the laminar rigid body rotation. The measuring 
point is in the developing region of the jet flow, 3 di- 
ameters away from the nozzle exit and is slightly off the 
centerline of the jet flow. Figure 4 shows the measured 
histogram P(Sv) of the velocity difference Sv(L,t) in the 
jet flow, when the separation L is fixed at L=0.5 mm (cir- 
cles) and L=0.8 mm (squares), respectively. It is seen 
that the measured P(Sv) has a dominant peak and its 
position changes with L. Because 5v(L,t) increases with 
L, the peak position moves to the right for the larger 



value of L. The solid curve in Fig. 4 is a Gaussian fit 
to the data points with L=0.5 mm. The obtained mean 
value of Sv(L,t) is (Sv) — 1.87 cm/s and the standard 
deviation a = 0.171 cm/s. At L=0.8 mm, the measured 
P(5v) peaks at the value (Sv) — 2.73 cm/s. 
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FIG. 4. Measured histogram P(8v) of the velocity differ- 
ence Sv(L, i) in the jet flow. The values of L are: L=0.5 mm 
(circles) and 0.8 mm (squares). The solid curve is a Gaussian 
fit to the circles. 

In the above discussion, we have assumed that each 
input fiber sees only one particle at a given time and 
the beat signal comes from two moving particles sepa- 
rated by a distance L. In fact, when the seeding density 
is high, each input fiber may see more than one parti- 
cle at a given time. The scattered light from these par- 
ticles can also beat and generate a self-beat frequency 
proportional to 5v(£q), where £q ~ 0.15 mm is the laser 
spot size viewed by each input fiber. ^ The self-beating 
gives rise to a small peak on the left side of the measured 
P(Sv). Note that the peak position is independent of 
L, because 5v(£q) is determined only by £q, which is the 
same for both measurements. It is seen from Fig. 4 that 
the cross-beating is dominant over the self-beating under 
the current experimental condition. 

The intensity correlation function g(r) is also used to 
analyze the beat signal. In the experiment, we measure 
the histogram P(5v) and g(r) simultaneously, so that 
Eq. (H) can be examined in details. Figure 5 shows the 
measured g(r) — 1 (circles) as a function of delay time r 
at L=0.5 mm. The squares are the self-beat correlation 
function G s (r) obtained when one of the input fibers is 
blocked. As shown in Fig. 4, the measured P(Sv) has a 
Gaussian form and thus the integration in Eq. (|4|) can 
be carried out. The final form of g{r) becomes 



g(r) = l + G s (r) 



i[q(Sv)n 



,-(qvr) 2 /2 



(5) 



The solid curve in Fig. 5 is a plot of Eq. (|5j) with b s = 0.5 
and b c — 0.13. The values of (Sv) and a used in the plot 
are obtained from the Gaussian fit shown in Fig. 4. It is 
seen that the calculation is in good agreement with the 
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measured g(r). 

The fitted value b s = 0.5 agrees with the expected 
value at I\ = 12- The value of b c would be 0.5 if the 
collected signals from the two input fibers were fully co- 
herent and the fiber-optic coupler mixed them perfectly. 
The fact that the fitted value of b c is smaller than 0.5 
indicates that the collected signals are not fully corre- 
lated. This is caused partially by the fact that in the 
present experiment the scattered light suffers relatively 
large number fluctuations resulting from a changing num- 
ber of particles in the scattering volume. These number 
fluctuations produce incoherent amplitude fluctuations 
to the scattered light and thus introduce an extra (addi- 



tive) decay to g(r) 
(proportional to the 



16[ Because the beam crossing time 
Deam diameter) is much longer than 
the Doppler beat time 1 / Au>2 (proportional to the wave- 
length of the scattered light), the slow decay due to the 
number fluctuations can be readily identified in the mea- 
sured g(r). This decay has an amplitude 0.4 and has 
been subtracted out from the measured g(r) shown in 
Fig. 5. 
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FIG. 5. Measured intensity correlation function g(r) — 1 as 
a function of delay time r at L=0.5 mm (open circles). The 
squares are obtained when one of the input fibers is blocked. 
The solid curve is a plot of Eq. (5). 

It is shown in Eq. (|^) that to accurately measure 
the mean velocity difference {Sv}, the beat frequency 
Aui 2 = q(Sv) must be larger than the decay rate T s ~ 
q5v(£o) for G s (t) and also larger than the decay rate 
r c ~ qa resulting from the fluctuations of the velocity 
difference. From the measurements shown in Figs. 4 and 



5, we find r s 



3.8 x 10 4 



which are indeed smaller than the beat frequency Acj2 — 
4.15 x 10 5 8 . Because g(r) contains a product of G s (t) 
and exp[— (qar) 2 / 2} [see Eq. (||)], its decay is determined 
by the faster decaying function. It is seen from Fig. 5 
that the decay of g{r) is controlled by G s (r), which de- 
cays faster than exp[— (qar) 2 /2]. It should be noted that 
in the measurements shown in Fig. 4, the beat signals 
are analogue ones and we have used a band-pass filter 
together with a LDV signal analyzer to resolve the beat 



frequency. Consequently, many low-frequency self-beat 
signals are filtered out. This low- frequency cut-off is ap- 
parent in Fig. 4. The measurements of g(r), on the other 
hand, are carried out in the photon counting mode, and 
therefore the measured g(r) is sensitive to all the self- 
beat signals as well as the cross beat signals. With a 
simple counting of particle pairs, we find that the proba- 
bility for cross beating is only twice larger than that for 
the self-beating. 



B. Local velocity measurements 

We now discuss the local velocity measurements using 
the new optical arrangement shown in Fig. 2(a). The ve- 
locity measurements are conducted on a freely suspended 
flowing soap film driven by gravity. Details about the ap- 
paratus has been described elsewhere, []T^]l9f and here 
we mention only some key points. 2% solution of deter- 
gent and water is introduced at a constant rate between 
two long vertical nylon wires, which are held apart by 
small hooks. The width of the channel (i.e., the distance 
between the two nylon wires) is 6.2 cm over a distance of 
120 cm. The measuring point is midway between the ver- 
tical wires. The soap solution is fed, through a valve, onto 
an apex at the top of the channel. The film speed v, rang- 
ing from 0.5 to 3 m/s, can be adjusted using the valve. 
The soap film is approximately 2-6 /im in thickness and 
is seeded with micron-sized latex particles, which scatter 
light from a collimated laser beam. The light source is an 
argon-ion laser having a total power of 1W. The incident 
laser beam is oriented perpendicular to the soap film and 
the scattered light is collected in the forward direction. 
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FIG. 6. Measured intensity autocorrelation function A(t) 
as a function of delay time r with the measuring time T=30 
ms. The inset (a) shows an enlarged portion of A(t) for small 
values of r up to r = 20 /is. The inset (b) shows the frequency 
power spectrum P(f) of the measured A(t). 

To measure the rapidly changing beat frequency Acl^, 
we build a fast digital correlator board for PC. |2(| 
With a fast sampling rate f s , the plug- in correlator 
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board records the time-varying intensity I(t) (number 
of TTL pulses from the photomultiplier tube per sam- 
ple time) over a short period of time T and then calcu- 
lates the (unnormalized) intensity autocorrelation func- 
tion, A(t) = (I(t + T)I(t)). Figure 6 shows an example 
of the measured A(t) as a function of delay time r with 
T = 30 ms and f s — 14.32 MHz. Because the burst sig- 
nal I(t) is a periodic function of t, the measured A(t) 
becomes an oscillatory function of r. The frequency of 
the oscillation apparent in the inset (a) is the beat fre- 
quency Au>i. The amplitude of the oscillation decays at 
large r. The inset (b) shows the power spectrum P(f) of 
the measured A(t); it reveals a dominant peak at 755.1 
kHz. The power spectrum is obtained using a fast Fourier 
transform (FFT) program. j2l| 
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FIG. 7. Measured intensity autocorrelation function A(t) 
as a function of delay time r with the measuring time T=50 
[is. The inset (a) shows the frequency power spectrum P(f) 
obtained by FFT. The inset (b) shows the frequency spectrum 
Q(f) obtained by the Scargle-Lomb method. 

To increase the sampling rate of the velocity measure- 
ments, one needs to keep the measuring time T for each 
A(t) as short as possible. The signal-to-noise ratio for 
A(t) decreases with shorter measuring time T and with 
lower mean photon count rate, which was ~ 1 MHz in the 
present experiment. It is found that the shortest useful 
measuring time T c is roughly 50 [is. For this value of T, 
A(t) is quite noisy and the corresponding peak in P(f) 
becomes less pronounced [see Fig. 7 and inset (a)]. It 
is worth mentioning that if one only wants to know the 
periodicity of a function, rather than its actual power at 
different frequencies, the Scargle-Lomb method [E2| is a 
better alternative to FFT. This method, which does not 
require evenly spaced sampling, compares the measured 
data with known periodic signals using the least-square 
fitting procedure and determines the relevant frequencies 
by the goodness of the fit Q(f). It can even utilize the un- 
even sampling to further increase the Nyquist frequency. 
As shown in Fig. 7(b), the Scargle-Lomb method can 
still clearly identify the periodicity of the signal, even 
when the power spectrum P(f) [Fig. 7(a)] becomes less 



reliable. The total time required for the measurement of 
the characteristic frequency is less than 1 ms. Using the 
correlator board together with an average speed PC (300 
MHz), we are able to conduct accurate measurements of 
the local velocity with a sampling rate up to 1 kHz. 



IV. SUMMARY 



We have developed new optical arrangements with two 
single-mode input fibers and a fiber-optic coupler to mea- 
sure the local velocity v(x) and the velocity difference, 
5v(£) = v(x + £) — v(x), between two spatial points sepa- 
rated by a distance I. The fibers and the coupler are po- 
larization preserving to guarantee a high signal-to-noise 
ratio. To measure the velocity difference 6v(£), the two 
input fibers are used to collect the scattered light with 
the same momentum transfer vector q but from two spa- 
tially separated regions in a flow. These regions are il- 
luminated by a single coherent laser beam, so that the 
collected signals interfere when combined via the fiber- 
optic coupler. The resultant light received by a pho- 
tomultiplier tube therefore contains the beat frequency 
Aui2 = q • Sv(£). We analyzed the beat signals using two 
different devices and compared the results. First, the in- 
tensity auto-correlation function g[r) was measured us- 
ing a digital correlator. Secondly, a standard LDV signal 
processor was used to determine the instantaneous beat 
frequency Aw 2 . With this device, 6v(£,t) can be ob- 
tained in real time. The technique can be further devel- 
oped to measure one component of the local flow vorticity 
vector d?(x, t)=Vx v(x, t). Q 

To measure the instantaneous local velocity itself, one 
needs only to reorient the two fibers so that they point to 
the same scattering volume. With this optical arrange- 
ment, we have three alternatives to measure a velocity 
component. They employ (i) an analog photodetector 
and a standard LDV signal processor (burst detector), 
(ii) a commercial photon correlator, such as that made 
by ALV, and finally (hi) a home-made digital correlator. 
This latter device completes a velocity measurement in 
less than 1 ms and is orders of magnitude cheaper than 
the other two alternatives. The new fiber-optic method 
has several advantages over the standard LDV and can 
be used widely in the general area of fluid dynamics. Be- 
cause only one laser beam is needed to obtain two veloc- 
ity components, a compact single-frequency diode laser 
can replace a large multi-frequency argon-ion laser. By 
eliminating the color and beam separation units in the 
standard LDV, the one-beam scheme is less costly to im- 
plement. With more optical fiber pairs and couplers, one 
can carry out multi-point and multi-component velocity 
measurements in various turbulent flows. 
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